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Abstract 



Q\ We investigate the entanglement dynamics of two initially entangled qubits interacting independently with two uncorrelated reser- 
voirs beyond the Markovian approximation. Quite different from the Markovian reservoirs [C. E. Lopez et al., Phys. Rev. Lett. 



(N 



101 (2008) 080503], we find that entanglement sudden birth (ESB) of the two reservoirs occurs without certain symmetry with 



respect to the entanglement sudden death (ESD) of the two qubits. A phenomenological interpretation of entanglement revival is 
Q^also given. 
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Introduction 

■y ' The entanglement dynamics of open quantum systems has 
attracted considerable interests over recent years lEl . By inves- 
^ tigating a system of two qubits interacting independently with 
O" 1 their corresponding vacuum reservoirs, Ref. yfl has pointed 
out the sudden termination of the entanglement initially owned 
by the two qubits at a finite interval, which is called entangle- 
;> ment sudden death (ESD). This intriguing phenomenon has re- 
cently been demonstrated experimentally by linear optics sys- 
£S terns JUS] and by atomic ensembles Jl. A great deal of fheo- 
■^j- retical inv estigati ons of ESD have so far been reported | 



i M El El Til El El El 



on the same model as in Ref. 



IIEIE 

'. m, it \ 



lUl . For example, based 
was found that ESD re- 



garding the qubits would always occur under the thermal and 
c^) squeezing reservoirs |S@]. Meanwhile, studies of ESD have 

been extended to multi -particle systems 1 1011 - which indicated 
• 1— 1 the delay of ESD time with more qubits involved. In addition, 

some studies have also showed that entanglement of qubits will 



revive in the case of a commonly shared reservoir!^ 



of independent non-Markovian reservoirs [13, 14, 15]. 



lilUJ or 



A question naturally arises: where has the lost entanglement 
gone when ESD occurs? An answer to the question has been 
given recently by |Q]. Considering the systems and the reser- 
voirs as a whole, Ref. 10] presented the lost entanglement to 
be transferred to reservoir degrees of freedom, which is called 
"entanglement sudden birth" (ESB) of the reservoirs. Intrigu- 
ingly, ESD of the systems and ESB of the reservoirs are of some 



symmetry and ESB might be occurring before, simultaneously 
with or even after ESD, depending on different initial states. 
However, the discussions in J7[] are restricted to the weak cou- 
pling regime under the Markovian approximation, which only 
works when the reservoir correlation time is small compared to 
the relaxation time of the system [20]. If the qubit is strongly 
coupled to its environment or the characteristic time of the sys- 
tem is shorter than the environmental correlation time 1 20.12 ill, 
such as in recent experiments with cavity QED system B22I1 or 
solid-state systems 123112411 . the situation would be much com- 
plicated and we have to give up the Markovian approximation 
in treatment. 

The present paper is focused on the study of entanglement 
dynamics without the Markovian approximation for two qubits 
interacting, respectively, with two independent reservoirs. For 
the initial condition with the two qubits entangled, we will study 
ESD and ESB, and show their different variation from the Marko- 
vian treatment. 

2. Physical model 

We consider a two-qubit system coupling to two uncorre- 
lated vacuum reservoirs at zero temperature. Since there is no 
interaction between the two pairs of qubit-reservoir, the dynam- 
ics of the whole system can be obtained simply from the evo- 
lution of the individual pairs Q7|, [l3[]. The Hamiltonian of the 
individual qubit-reservoir pair under the rotating wave approx- 
imation is given by (ft = 1), 
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H = <y o" + cr_ + 2^ (±> k a\ ak + 2^ 8k(cr-al + cr+a k ), (1) 



k=\ 



k=\ 



where loq is the resonant transition frequency of the qubit be- 
tween levels |0> and |1>, with cr + = 1 1 > <0| and cr_ = |0><1|, 
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<jL>k, a + k (a.k), gk are the frequency, creation (annihilation) oper- 
ator and the coupling constant for the A:th mode of the reser- 
voir. For simplicity in our treatment, only one excitation of 
the reservoir will be considered. Consequently, the initial state 
|Y(0)> = C (0)|l>„|0> r + Eli C*(0)|0>Jl*> r will evolve to 



Using Eq. (5), we could straightforwardly reach the reduced 
density matrices for the two qubits, 



p, 192 (t) = 
-dco 



\y(t» = c «|i> e |o> r + J] c t (t)|o> 9 |iit> r , 



(2) for the two reservoirs, 



k=] 



where |0),- = 



|0jt) r , and |ljfc) r , an abbreviation of ( 



h N 



) k=l \-',,/f I .„/,'. c,,.>.„^ , /= | / = A 

is the reservoir state with one excitation in the Mi mode and 
other states in vacuum. Assuming the initial conditions Co(0) = 
1 and Ck{0) = 0, we could obtain the following equation by 
solving the Schrodinger equation of Eq. (2), 



|0,->,|U>r, 



Pr 1 r 2 (') = 



C (f) 



r 

-/ 



dt Fit - t )C (f ), 



(3) 



where the correlation function F(t — t ) in the limit of N — > 
oo is of the form F(t - t) = j da>J((j)e Ka °- b ' )(/ -t') with J(a>) 
the spectral density of the reservoir. In what follows, we will 
consider in our calculation the Lorentzian spectral distribution 



J(oj) = 



— 12011 . in which yo — T o an( l J — T r , with 



2k (a, -c 1 )) 2 +y 2 i±aai>-"» /u ' 

to and r, the qubit relaxation time and the reservoir correlation 
time, respectively H20TI . In a non-Markovian regime, yo > y/2, 
i.e., the reservoir correlation time is longer than the relaxation 
time of the qubit, and yo < y/2 means a Markovian regime 
[20]. As a result, in the non-Markovian regime, the probability 
amplitude can be easily solved as 



C (f) = e 



COS(y)+ -Sin(y) 



(4) 



with T = Vr(2yo - y). If we set C(f) = y/l - C (t) 2 , Eq. (2) 
can be rewritten as 

|Y(0>=Co(OH>,|8>r + C(OIO>jI>r, (5) 

with |l> r = [2f =1 C k (t)\h) r ] /C(t) (N oo). 

3. Entanglement dynamics 

We will employ concurrence B25I1 to assess entanglement. 
In the case of the density matrix with the "X" type 





a 



w" 



z 




w 



d 



(6) 



the concurrence could be easily calculated by [2 

C (p) = 2 max{0, \z\ - <^d, |w| - V^}. (7) 

We consider that the two qubits are initially entangled and 
the two reservoirs are in zero-temperature vacuum states, which 
is described as 



P(0)=- 





1 
1 









1 -a 



>|0> ri |0> r 



(8) 



aCjj(!)C z (t) + C z </> 

c 2 (0 



C*(0 



aC 4 (() + 2C 2 (0+l- 



(9) 



ad(t) + 2C 2 (0+l- 



(10) 



for the qubit- 1 and reservoir- 1 , 



l + «)C^(f) 
l-»)C (t)C(/) 



«)C «)C(0 

+ a)C 2 (t) 



(ID 



and for the qubit- 1 and reservoir-2, 



P3l^2 CO = 

oc2(r)C 2 (/) 



C 2 (l) + aC*(t) 
C (r)C(f) 



C (t|C(0 
C 2 (0 + oC^CO 



(12) 



As Eqs. (9)-(12) are all of the "X" types, the concurrence of 
the partitions q\-qi, r\ —r%, q\-r\ and q\ — T2 can be calculated 
directly using Eq. (7). In addition, Eq. (7) also implies that the 
ESD time of the qubits or the ESB time of the reservoirs can be 
obtained from \z\ - yfad = and \w\ - yfbc = 0. However, as we 
could not solve them in an explicitly analytical way, numerical 
calculations are given below instead. 

4. Numerical results and discussions 

Fig. 1 shows the entanglement dynamics of q\ - qi, r\ - r%, 
q\ — r\, andgi —r% in the non-Markovian regime with y = O.lyo. 
In Fig. 1(a), there are three regimes for entanglement variance 
regarding the two qubits: for a < 0.3, no ESD will occur; for 
0-3 ^ a < 0.4 there will be entanglement revival after ESD 
takes place; as for a > 0.4 entanglement will vanish forever at 
a finite interval lfl3ll . In contrast to the behavior regarding the 
two qubits, there are only two a-dependent regimes for entan- 
glement variance regarding the two reservoirs: ESB or no ESB, 
as shown in Fig. 1 (b). But no matter in which regime, the 
two reservoirs will finally be entangled in a time-dependently 
oscillating manner. The entanglement transfer^ is done step by 
step from the qubits to the reservoirs, which is reflected in Figs. 
1(c) and 1(d) with the entanglement variance regarding the par- 
titions q\ — r\, q\ — r% etc. . The partitions q\ - r\ , q\ - behave 
as something like relays, which get and release entanglement 
in an oscillating way. Meanwhile the r\ — entanglement is 
enhanced also in the oscillating way. 



2 To avoid misunderstanding, we would like to mention that the term "en- 
tanglement transfer" used in the present paper means the appearance of a lost 
entanglement in some other degrees of freedom, e.g., qubit-reservoir, reservoir- 
reservoir etc. due to complicated couplings between the qubits and the environ- 
ment. This is quite different from in Refs. 1 27, 28, 29, 30, 31, 32] etc., where 
"entanglement transfer" is referred to the movement of entanglement from the 
flying qubits to static qubits by some operations. 




Figure 1: (Color online) Concurrence as a function of a and yot for Eqs. (9)-(12) in the non-Markovian regime with y = O.lyo- 




Figure 2: (Color online) Concurrence as a function of yot with a = 0.35 for Eqs. (9)-(12) in different coupling intensities where the green (grey) solid lines mean 
the Markovian case with y = 5yo; the blue curves correspond to the non-Markovian cases with y = 0.1 yo (solid lines) and y = 0.05yo (dashed lines) respectively. 
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To see more clearly and also to make comparison with the 
Markovian treatment, we have plotted Fig. 2 for the entan- 
glement dynamics with a = 0.35 {t£SD > fee) - In the non- 
Markovian regime, e.g., y = 0.05yo, when ESD occurs at yot 
8, the entanglement of the qubits has been transferred to other 
partitions, e.g., r\ - T2, q\ — r\, q\ — r2 etc. . In more details, we 
may first visualize that within an interval yot e [0, 5.5), in which 
the entanglement transfer from qubit-qubit to qubit-reservoir is 
greater than that from qubit-reservoir to reservoir-reservoir de- 
grees of freedom, and the entanglement of qubit-reservoir and 
reservoir-reservoir are both increasing. Then from the critical 
point yot ^5.5, the qubit-reservoir entanglement begins to be 
decreasing until yot 11, at which the entanglement of qubit- 
reservoir has been totally transferred to the reservoir degrees of 
freedom. However, in contrast to the case of Markovian reser- 
voirs with the entanglement fully transferred to the reservoirs 
at one time, due to the memory effect of the non-Markovian 
reservoirs, entanglement will partially be recoiled to the qubit- 
reservoir degrees of freedom, or even back to the qubit-qubit 
degrees of freedom. As a result, revival of entanglement takes 
place HUHl- This back and forth process will repeat until the 
entanglement of the qubits have been completely moved to their 
corresponding reservoirs. Finally, the reservoirs are entangled. 
Moreover, Fig. 2 also presents that, the stronger the qubits cou- 
pling to their reservoirs, the more evident the non-Markovian 
memory effects. 

5. Conclusions 

In summary, we have investigated two-qubit entanglement 
affected by the decoherence from non-Markovian reservoirs. In 
contrast to the smooth decay of the qubit entanglement in the 
case of Markovian reservoirs, the entanglement of the qubits 
will revive after ESD occurs in the non-Markovian case. Mean- 
while, by means of the qubit-reservoir partitions as the relays, 
the reservoirs will be asymptotically entangled in an oscillat- 
ing manner, which is a good interpretation of entanglement re- 
vival of qubits in the non-Markovian reservoirs. We argue that 
our present study would be useful for quantum information sci- 
ence. For example, the solid-state qubits in systems with strong 
correlation always experience strong decoherence 112 lL I23L 12411 . 
So quantum information processing on any qubit in such sys- 
tems should be affected by detrimental influence from the non- 
Markovian reservoir. Moreover, fast logic gating is necessary 
in quantum information processing. As the characteristic time 
of the system during the fast gating is very short, we have to 
consider the non-Markovian effect from environment in assess- 
ing the operational fidelity. We have noticed recent studies on 
non-Markovian reservoir from the microscopic way 11331 13411 . 
Our present work from the phenomenological viewpoint, how- 
ever, could explain some experimental observations of entan- 
glement loss regarding qubits more straightforwardly. On the 
other hand, entanglement is the unique resource in quantum in- 
formation science. As no system could be completely isolated 
from external noise, it is of interest to understand what happens 
in the entanglement transfer from the system to the environ- 



ment, although whether the reservoirs are entangled or not is 
actually of no practical application with current technology. 
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